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1 Introduction 



The apparent freedom in choosing the renormahsation scale and scheme for per- 
turbative calculations of observables in quantum chromodynamics (QCD) introduces 
theoretical uncertainties which, if taken literally, prohibit absolute predictions be- 
yond a qualitative level. The renormahsation scheme dependence can be solved by 
using a fixed reference scheme or, equivalently, by relating measurements of different 
observables to each other. The MS scheme is the simplest choice from a calculational 
point of view but the question then arises if there exists a preferred scheme which 
is optimal from a physics point of view. A closely related question is how to choose 
the renormahsation scale which is important since most QCD observables are only 
known to next-to-leading order (NLO) where the renormahsation scale dependence is 
still sizable. 

Another problem with perturbative QCD predictions is that the series is in fact 
asymptotic, i.e. after a given order the higher order contributions start to increase 
and make the series divergent. The most prominent source for this asymptotic be- 
haviour is due to so-called renormalons which make the higher order coefficients grow 
factorially 

This talk presents an alternative approach which avoids, or at least minimizes, the 
problems outlined above by using conformal expansions and the closely related skele- 
ton expansion. The presentation is mainly based on which also contains a complete 
list of references. The relation between the skeleton expansion and the Banks-Zaks 
expansion as well as the BLM scale-setting method by Brodsky, Lepage, and 
Mackenzie and its generalizations is also discussed. 

2 Conformal relations 

For definiteness and simplicity the discussion will be limited to single-scale space- 
like observables in massless QCD, but the approach can also be generalised to time- 
like and multi-scale observables. The perturbative expansion for such a single-scale 
observable can be written as. 



where = —q^ is the (space-like) physical scale, /x^ is the renormahsation scale, 
and is the next-to-next-to-leading order coefficient in the renormahsation group 
equation for the coupling. 
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where a = as/n. 

The truncation of the perturbative expansion at order introduces a renormal- 
isation scale and scheme uncertainty of order a^"*"^. In addition the perturbative 
coefficients i?„ will asymptotically grow factorially due to renormalons, Rn ~ n!/5g . 
This should be contrasted with the situation in the conformal (scale-invariant) limit 
where da/dlog{fi'^) = 0. In this case there is no scale-ambiguity, and the coefficients 
Rn are free of factorial growth due to renormalons. The only remaining problem is 
the scheme uncertainty which can be circumvented by relating observables to each 
other instead of trying to make absolute predictions. 

Before continuing it is useful to recall the concept of an effective charge |^ which 
collects all perturbative corrections to an observable. An observable R{Q^) can then 
be written in terms of the effective charge CLRiQ"^) as, 

R{Q') = i?QPM(g') + R<,{Q')aR{Q') 

where 

and the perturbative coefficients = Ri/Rq. 

The most celebrated example of a conformal relation between observables is the 
Crewther relation |jTO|,|lT|,|^ between the Adler D-function (a^) and the polarized 
Bjorken sum-rule for deep inelastic scattering (a^J, 

(1 + a^)(l - OgJ = 1 . 

Thus, the Crewther relation is simply a geometric series to all orders and there is no 
growth of higher order coefficients. The effective charges ao and a^^, which appear 
in the relation, are defined by. 



Dm = Q'^^ = iVcEe? [l + aniQ' 







/ 

' gfix,Q'')-g'lix,Q^) dx = I - a^^iQ' 



Qgv 



where Il{Q^) is the hadronic correction to the vacuum polarisation of the photon, the 
spacelike continuation of Re+e-{s). 

In general, conformal relations between two arbitrary observables A and B can be 
written as, 

n 

where, as is evident, the conformal coefficients c^^ depend on which two observables 
that are related. Of course, in real life the QCD coupling is scale-dependent. Even 
so, the notion of conformal coefficients is still useful as will be shown below. The 
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main advantage is that by identifying the conformal part of the ordinary perturba- 
tive coefficients it is possible to treat all the running coupling effects separately and 
thus keeping the coefficients free from factorial growth due to renormalons which are 
instead resummed in the running of the coupling. 



3 The skeleton expansion 

The skeleton expansion |jl2| organizes the perturbative series in terms of contribu- 
tions to fundamental skeleton graphs. A skeleton graph is defined by the requirement 
that the fundamental vertices and propagators contain no substructure. One example 
of an ordinary Feynman diagram and the corresponding skeleton graph is shown in 

Fig.g. 




Figure 1: Example of an ordinary Feynman diagram (left) and the corresponding skeleton 
graph (right) in QED. 



In QED, the skeleton expansion is straight-forward to construct thanks to the 
basic Ward identity, Zx = Z2, from which it follows that charge renormalisation is 
given by photon propagator renormalisation (Z^). The coupling a that appears in 
the skeleton expansion is the Gell-Mann Low coupling which resums the Dyson series 
of the one-particle irreducible photon self-energy 11, 



1 - n(Q2) 



The radiative corrections to the one-photon exchange skeleton graph, such as illus- 
trated in Fig. 1^, can then be written as an integral over the running coupling. 




where 0o is a momentum distribution function which has been normalised to 1 for 
convenience. (In the above example 0o (^) = ^{^'^ ~ Q"^)-) 
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Adding the contributions from one-, two-, three-photon exchange etc., an effective 
charge can be written as 



„ / h.2 h.2 U2\ JU2 JU2 JU2 



"-1 "-2 



where (pi are the momentum distribution functions (normahsed to 1) and the q are 
the conformal coefficients in the skeleton scheme. For simphcity the above expression 
has been written including just one skeleton at each order but in general there can be 
several different skeletons which contribute at the same order. For comparison, the 
conformal theory gives aji{Q'^) = a + CiO? + 020?' + ■ ■ ■. 

Another important property of the skeleton expansion is that each term in the 
expansion is renormalisation scheme and scale-invariant by itself. In addition the 
skeleton coupling is gauge-invariant. The skeleton expansion thus provides an alter- 
native way of writing the perturbative series for an observable in which each term 
is given by one or several integrals over the running coupling. One complication 
of the skeleton expansion is that in general one needs a diagrammatic construction 
to identify the different skeletons. However, at low orders this requirement can be 
bypassed. 

In QCD, the existence of an all-order skeleton expansion has so far not been 
proved. The basic complication arises from the gluon self-interactions and the related 
difference between gluon-propagator and charge renormalisation. Nevertheless it is 
reasonable to assume that something similar to the skeleton expansion in QED can 
also be constructed for QCD. In fact, the so called pinch technique [l^ provides a 



realisation of the skeleton expansion in QCD at the one-loop level. As an example 
Fig. H illustrates how the three-gluon vertex is divided into a pinch part which con- 
tributes to the renormalisation of the effective propagator and a non-pinch part which 
contributes to renormalisation of the "external" vertex. 



Figure 2: Illustration of the subdivision of the three-gluon vertex into a pinch part and a 
non-pinch part using the pinch technique 



4 



In this way the pinch technique arrives at the following QED-like Ward identities, 



1+ 



- 1 :rCi 

e 4 

—Ca- -TpNp 
12 3 



1 + -(3oa 

e 



such that all the one-loop running coupling effects are contained in the effective gluon 
propagator. The coupling defined by the pinch technique has a simple relation to the 
MS scheme, 



a{Q^) = aMs(/i^) + 



+1 



Recently there has been progress in extending the pinch-technique to two loops |]T1| 
and this may eventually lead to an extension of the skeleton expansion in QCD to 
two loops as well. Another possibility may be to use light-front quantization of QCD 
in light-cone gauge [p!5[] . 



4 Identifying conformal coefficients 



Given the advantages of the skeleton expansion compared to the standard per- 
turbative expansion, it is instructive to consider the following simplified ansatz for 
QCD as a starting point for further investigations: assume there is only one skeleton 
coupling, that there is only one skeleton graph at each order in a, and that the depen- 
dence on the number of light flavours (Np) can be used to identify the non-conformal 
parts of the perturbative coefficients. Given these assumptions the first conformal 
coefficients in the skeleton expansion can be obtained from the perturbative ones in 
the following way 0. 

The starting point is the skeleton expansion of an effective charge given by Eq. (|l|). 
Next the skeleton couplings a(fc^) under the integration sign can be expanded in the 
coupling d{Q^) using the solution to the renormalisation group equation. 




Pi lo, 




+ /5o^log^ 




Inserting this into Eq. (|l|) then gives 
where 

log 



= a(g2) + (ci + /5o0[,')) a\Q^) + (c2 + cM? + Pi^'i^ + Po^^o^) a'(Q') + • 
J-""* are log-moments of the momentum distribution functions. 



k^ 



5 



This can now be directly compared with the standard perturbative expansion, 



which gives the relations 



ri = ci + /3o0o 



Based on the Np dependence of the perturbative coefficients fj it is thus possible 
to identify in a unique way the conformal coefficients Ci and C2 as well as the log- 
moments ^o^'*, 0o^\ and (This follows since the coefficients fj are polynomials 
in Np of order i.) In fact, given the assumptions made, it is possible to decompose 
the perturbative coefficients up to order without any additional information. At 
higher orders the Np dependence alone does not provide enough information even 
with the simplifying assumptions that have been made. 

In general there are several ways in which the assumed ansatz can break down. 
Most notably, at higher orders there are skeletons which are A'"i?-dependent by them- 
selves. In contrast to QED where the iVir-dependent skeletons (such as the light-by- 
light scattering diagrams) can be easily identified based on the dependence on the 
external charge there is in general no such simple identification possible in QCD. 
Another complication is that there may be more than one skeleton at each order. To 
resolve these two problems one will need an explicit diagrammatic construction of the 
skeleton expansion. It may also be the case that the skeleton expansion in QCD can 
only be systematically extended to all orders by having several skeleton couplings. 
However, even if some of the assumptions that have been made are wrong, it may still 
be true that the general properties of the ansatz are valid. This includes the prop- 
erty that running-coupling effects can be associated with different skeleton graphs in 
a renormalisation-group-invariant way, and that the skeleton coefficients are confor- 
mal. In practice there is usually no problem in identifying the skeleton structure at 
next-to-leading order but special care has to be taken as will be discussed below when 
the application of BLM scale-setting to the thrust-distribution in e"'"e~-annihilation 
is re-examined. 



5 Relation to the Banks-Zaks expansion 

As already realised at the time of the discovery of asymptotic freedom, perturba- 
tive QCD has an perturbative infrared fixed-point |T^ (A;^ — 0), 





din k"^ 
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in the so called conformal window 8 < Np < 16 since for this range of Np the first two 
terms in the /3-function have opposite signs, /3o = x ~ ^ ^ ~ ^ ~ ^ 

0. 

If the coupling at the fixed point Opp is small, such that perturbation theory is 
still applicable, then it can be written as a so called Banks-Zaks expansion in the 
parameter oq = -po/ Pi\i3o=o = w/^o, 

app = Oo + Vial + • • • , 

where the coefficients can be calculated from the higher order terms (/32 etc.) in 
the /3-function. 

In the same way an arbitrary effective charge can also be expanded in oq. 
Starting from the ordinary perturbative expansion the coefficients can be rewritten 
in terms of Oq using the polynomial A'^p-dependence, 

aniQ^) = a{Q^) + (ri,o + n^iao)a\Q^) + (r2,o + r2,iao + r2,2a^)a^(Q^) + • • • • 

From this it follows that it is also possible to get a relation between the fixed-point 
value of the effective charge and the coupling Opp. Taking the limit — > 
(assuming that this is well defined) and inserting oq = ctpp + wiOpp + ■ ■ ■ gives the 
fixed point relation, 

= Opp + ri^oflpp + (r2,o + ^^i.^ctpp H • 

Comparison with the conformal coefficients obtained from the skeleton decomposition 
of the perturbative coefficients shows that, if Cpp is identified with the skeleton cou- 
pling then, they are indeed the same, i.e. ri o = Ci and r2,o + '"i,! = C2 etc. Thus, the 
conformal coefficients in QCD can also obtained from the Banks-Zaks expansion by 
analytically continuing the number of light quark flavours into the conformal window 
and taking the infrared limit p[. 

6 Connection to BLM scale-setting 

Once the conformal coefficients have been identified one also has to evaluate the 
corresponding skeleton integrals. For the leading skeleton this can be done using 
the momentum distribution function calculated in the large /?o-approximation. At 
the same time the associated renormalon ambiguity indicates the form of the non- 
perturbative corrections in terms of power-corrections. The combination gives a 
framework for analysing the renormalon resummation and the non-perturbative cor- 
rections together [0,^. An alternative is to approximate the skeleton integrals by 
using BLM scale-setting |^,^ as will be discussed below. 
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The starting point is the skeleton expansion of the effective charge in question 
where each integral is evaluated using the mean value theorem (MVT) in the following 

way, 



+ c./a(.DaM0.(|,|4)f 

(MVT) = a{kl)+cia^{kl)+C2a^{kl) + --- 

The "BLM" scales ko, ki, k2, etc. are uniquely determined by requiring a one-to-one 
correspondence between the skeleton integrals and the terms in the "BLM" series [0]. 
In other words ko depends only on (j)Q, ki on 0i, and so on. Thus there is no ambiguity 
in determining the scales as is the case for commensurate scale relations Ex- 
panding the couplings d{k'^) in terms of a(Q^) under the integration sign the "BLM" 
scales are obtained as a perturbative series in the skeleton coupling with the coeffi- 
cients given in terms of the moments of the distribution functions, 




9^ - 

kn 







}f - ( 





/5oa(fco) + ■ ■ 



variance 
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It is important to realize that this provides a systematic improvement of the original 
BLM-scale, /cq^blm = exp ^o^j. In the lowest order approximation the scale fco 
is simply given by the mean of the momentum distribution as indicated above. By 
going to higher orders one then takes into account the variance of the distribution 
and so on. This corresponds to performing the skeleton integral with successively 
improved approximations to 0o- 

Given the conformal expansions of two observables in the skeleton scheme it is also 
possible to eliminate the skeleton scheme and get a direct relation between the two 
observables - a so called commensurate scale relation (GSR). From renormalisation 
group transitivity it follows that the coefficients in the commensurate scale relation 
are also conformal and thus free of factorial growth due to renormalons. However, 
there is no clear interpretation of the scales that appear in the GSRs, and in addition 
there is no unique scale setting procedure as has been already mentioned. 



7 Re-examining BLM scale-setting for thrust 

The new insights gained from the relation between the skeleton expansion and 
BLM scale-setting makes it interesting to re-examine BLM scale-setting for event 
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shape observables in e~^e annihilation [^. In the following the thrust distribution 



will be considered as a concrete example but general criteria for the applicability of 
BLM scale-setting will also be given. 

Thrust is an event shape observable defined by, 

EifiT- Pi 

1 = max ■ 



Ei \Pi 



where the sum runs over all particles in the final state. The thrust-axis ut is varied 
until the maximal value for T is obtained. An event with two narrow back-to-back 
jets corresponds to T = 1 whereas the minimal thrust value T = 0.5 is obtained for 
an event with isotropic distribution of particles as illustrated in Fig. ||. 

two-jet isotropic 



T=1.0 T=0.5 



Figure 3: The values of thrust for an event with two narrow back-to-back jets (left) and an 
event with isotropic distribution of particles (right). 



In the quark parton model the thrust distribution is a delta-function at T = 1. 
The leading order QCD-corrections have been calculated analytically |^ whereas the 
next-to-leading order QCD-corrections have only been calculated numerically [pl| ,p2 



From the definition of thrust one expects that at leading order there is only one skele- 
ton which contributes and that all the iVi?- dependence at next-to-leading order is from 
running coupling effects. Thus BLM scale-setting should be straight forward. There 
is however one possible complication, namely the non-inclusiveness of the definition. 
The easiest way to see this is that at leading order thrust can have values in the range 
2/3 < T < 1 whereas at next-to- leading order the range is given by 1/^3 < T < 1. 
Thus, if the next-to-leading order A^j?-dependence is non-zero for l/^/3 < T < 2/3 
then this cannot be attributed to the leading skeleton. However, as will be shown 
below the problems for the case of thrust are minimal. 

At next-to-leading order the BLM series for the thrust-distribution can be written 

as, 

-^^(^,T) = 5(l-T) + i?o(T)a(A:oVM(^,T))+i?i(T)a2(fco%,^(.,T)) , (2) 

where the skeleton coupling has been identified with the pinch technique coupling, 
-Ri(T) is the conformal coefficient in the pinch scheme, and the BLM-scale fcosLM is 
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Figure 4: (a) The next-to-leading order conformal coefficient compared to tlie standard 
MS coefficient RiusifJ-'^ = s,Np = 5) and the leading order (scheme-invariant) coefficient 
Rq. (b) The BLM-scale for ^/s = Mz- For both figures the points show the numerical values 
that have been calculated and the lines are fits to these points taking into account the known 
logarithmic terms. In (b) the dotted line show the approximation A;o,blm — 1.4(1 — T)^/s. 



used to approximate the unknown scale which should appear in the ^i(T)-term. It 
is important to realize that the BLM-scale /cq blmI"^' ^) ^ function of both kinematic 
variables, s and T. In addition the BLM scale is undefined for T < 2/3 where Rq 
vanishes. 

The expansion given above should be compared with the standard MS expansion 
using fi'^ = s, 

1 d(T^ 

--— (s, T) = 5(1 - T) + i?o(T)aMs(s) + RiMf^' = s, Np, T)aUs) • 
a al 

The leading order coefficient Rq is scheme-invariant and thus the same in both expan- 
sions. However, the next-to-leading order coefficient Ri is very different in the two 
cases as is illustrated in Fig. ^ which shows the conformal coefficient compared to 
the standard MS coefficient i?i^Ms(/^^ = -s, Np = 5) and the leading order coefficient 



Rq. The coefficients have been calculated numerically using the Beowulf program [22 



which is shown as points in the figures. The lines are fits to this points taking into 



account the know logarithmic parts of the coefficients [25]. 

From the figure it is clear that the next-to-leading order coefficient is large com- 
pared to the leading order one in both cases. However, the conformal coefficient is 
more stable over a large range of T (when multiplied with (1 — T)) except for T — >■ 1 
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Figure 5: (a) Fit to OPAL data using the fixed order BLM series. The fuh Une corresponds 
to the range fitted (0.7 < T < 0.95). (b) Value of running couphng as extracted from 
OPAL data in the range 0.7 < T < 0.95 at the corresponding BLM-scale. 

where it becomes negative. This is the Sudakov region which can only be properly 
treated by resumming all singular terms in the Sudakov form-factor. Another im- 
portant feature which is clear from the figure is that the non-conformal part of Ri 
more or less vanishes for T < 2/3, which is a good indication that the Np dependence 
can indeed be used to separate the conformal and non-conformal parts and that the 
problems with non-inclusiveness are only minor (see also p3| , p!8| , p^ ) . This property 
is different for other event shape observables depending on how they are defined. For 
example, oblateness is defined as the difference between an observable that starts at 
order as and one that starts at order a^. As a consequence there are Np depen- 
dent contributions to the next-to-leading term which do not come from the leading 
skeleton. This could also explain why BLM scale-setting seems to fail for some event 
shape observables [0. 

Fig. H also shows the resulting BLM-scale /co.blmI-s, 7") for the case i/i = Mz- 
From the figure it is clear that the scale vanishes as T — > 1 which is reasonable since 
the available phase-space for gluon emission vanishes in this limit. For comparison 
the figure also shows the approximation ^cblm — 1.4(1 — T)y/s which gives an overall 
good description of the T-dependence. The scale can also be understood physically 
as the transverse momentum which approximately scales as (1 — T)y/s for a three-jet 
configuration with one of the jets being much less energetic than the other two, i.e. in 
the soft limit. For T — >^ 2/3 the BLM-scale grows rapidly since the i?o — ^ but even 
at T = 0.69 (the point with the smallest T-value shown in the figure) the BLM-scale 
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is still smaller than ^/s which should be true in a physical scheme following from the 
mean value theorem. 

For illustration, the fixed order BLM expression for the thrust distribution given 
by Eq. has been fitted to data from the OPAL collaboration [^| at ^/s = Mz 
in the range 0.70 < T < 0.95 using a two-loop running coupling. The result of 
the fit, which is shown together with the data in Fig. ^, corresponds to the value 
aMs(M|) = 0.117. (To translate the fit into a value for aMs(M|) the commensurate 
scale relation, a(e^/'^M|) = aMs(-^|)+o|jg(M|), was used.) For comparison, using the 
fixed order MS expression gives aMs(M|) = 0.143. This illustrates the importance of 
taking running coupling effects into account. However, it should be kept in mind that 
a complete analysis should also include the Sudakov form-factor and non-perturbative 
effects. 

It is also possible to see the running of the coupling function of the BLM- 

scale A;o,BLM directly from the data. For each data point Eq. (|D is a simple second 
order equation which can be solved for d = Os/tt. The resulting values of obtained 
in this way are shown in Fig. ^ as a function of the corresponding BLM-scales. (The 
figure only shows the points that where used in the fit. For larger values of T the 
next-to-leading order coefficient -Ri is negative and for smaller values of T the next- 
to- leading order correction is larger than 100%.) Thus, even though the experiment 
is done at a fixed energy, it is still possible to observe the running of the coupling. 



8 Conclusions 

The standard perturbative expansion of observables in QCD is plagued by renor- 
malisation scheme and scale ambiguities as well as higher order coefficients which grow 
factorially due to renormalons. In this talk I have presented an alternative approach 
which avoids, or at least minimizes, these problems by using conformal expansions, 
especially the skeleton expansion. 

In contrast to the ordinary perturbative expansion the skeleton expansion is free of 
renormalisation scheme and scale ambiguities and the coefficients are free of factorial 
growth due to renormalons. Presently the pinch technique provides a realization of 
the skeleton expansion in QCD at next-to-leading order but it is not known whether 
an all-order expansion exists or not. Even so, the skeleton expansion has important 
phenomenological consequences. 

The leading skeleton integral makes it possible to include non-perturbative effects 
in a consistent way which takes into account the arbitrariness of the definition of 
perturbation theory. The renormalon ambiguities which appear in the evaluation 
of the leading skeleton integral can be used to parametrize the non-perturbative 
contributions in the form of power-corrections. 

By making a simple ansatz for the skeleton expansion in QCD the first steps in 
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making a more systematic study of it properties have been taken 0. One result of 
this study is that the conformal coefficients coincide with the ones obtained in case 
QCD has a perturbative infrared fixed-point (the Banks-Zaks expansion). 

The skeleton integrals which appear in the skeleton expansion can also be ap- 
proximated by the BLM-scale setting method and its generalisations. Requiring a 
one-to-one correspondence between the BLM-scales and the skeleton integrals gives 
a unique prescription for setting the scales in contrast to the situation for com- 
mensurate scale relations. The connection between the skeleton expansion and BLM 
scale-setting also gives new criteria for the applicability of the latter. 
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